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ABSTRACT 

Control -configured vehicle technology has Increased the demand for 
detailed analysis of dynamic stability and control, handling and ride 
qualities, and control system dynamics at early stages of preliminary 
design. An approximate, but reasonably accurate, set of ejuatlons of 
motion are needed for these early analyses. Such a formulation Is 
developed for the longitudinal dynamics of elastic airplanes. It makes 
use of only rigid-body aerodynamic stability derivatives In formulating 
the forces and moments due to elastic motion. Verification of accuracy 
using data for the B-1 airplane shows very good agreement. Frequencies 
and damping ratios uf the coupled modes corresponding to complex roots of 
the characteristic equations agree closely with four symmetric elastic 
modes Included. 
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WOMLNCLATURE 
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■ mean aerodynamic chord (ft) 

■ drag coefficient 

• lift curve slope stability derivative 

■ lift coefficient due to elevator deflection stability deriva- 
tive. 

» pitching moment coefficient due to angle of attack stability 
derivative 

■ pitching moment coefficient due to wing dowfiwash on tall 
stabll Ity derivative 

■ pitching moment coefficient due to pitch rate stability 
derivative 

■ pitching moment coefficient due to elevouwi’ deflection stabil- 
ity derivative 

■ 1th e'lastic mode aerodynamic force coefficient in y-directlon 
due to plunge velocity of C.G. (1/sec) 

■ 1th elastic mode aerodynamic force coefficient in z-directicn 
due to jth mode generalized displacement (l/sec‘) 

» ith elastic mode aerodynamic force coefficient In z-direction 
due to jth mode generalized velocity (1/sec) 

= 1th elastic mode aerodynamic force coefficient in z-directlon 
due to elevator deflection (ft/sec^) 

“ ith elastic mode aerodynamic force coefficient In z-directlon 
due to vertical gust velocity (1/sec) 

“ mass moment of inertia about y-axis (slug-ft^) 

* total airplane mass (slugs) 
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M 


w 




M. 


M. 




M" 


M. 


qg(t) 

S or S, 


Uo 

w(x,y,t) 

Wg(t) 


w 


aerodynaric pitching moment stability derivative due to plunge 
velocity of C.G. (rad/ft-sec) 

aerodynamic pitching moment stability derivative duf to down- 
wash from wing to tail (rad/ft) 

aerodynamic pitching moment stability derivative due to pitch 
rate (1/sec) 

aerodynamic pitching moment stability derivative due to 
elevator deflection (l/seo^) 

aerodynamic pitching moment coefficient due to ith elastic 

mode generalized displacement (rad/ft-sec^) 

aerodynamic pitching moment coefficient due to ith elastic 

mode generalized velocity (rad/ft-sec) 

aerodynamic pitching moment coefficient due to ith elastic 

mode generalized acceleration (rad/ft) 

j I m(x,y)(>^^(x,y)dxdy 1th elastic mode generalized mass 

^ ^ (slugs) 

pitch gust velocity (rad/sec) 

wing planform reference area (ft^) 

horizontal tail planform reference area (ft^) 

trim flight velocity (ft/sec) 

local plunge velocity in z-direction (ft/sec) 

vertical gust velocity at C.G. in negative z-direction (ft/sec) 

aerodynamic force stability derivative in z-direction due to 

plunge velocity of C.G. (1/sec) 

aerodynamic force stability derivative in z-direction due to 
elevator deflection (f t/rad-sec^) 
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C^(t) 

e(x,y,t) 

Po 

u»1 

♦^(x.y) 

>j(x,y) 


aerodynamic force coefficient In z-directlon due to 1th 
elastic mode generalized displaceiiient (1/sec^) 
aerodynamic force coefficient In z-d1rect1on due to 1th 
elastic mode generalized velocity (1/sec) 
local angle of attack (rad) 
elevator deflection (rad) 

1th elastic mode structural damping ratio 

1th elastic mode generalized displacement In z-directlon (ft) 

local pitch angle (rad) 

free stream air density (slugs/ft^) 

free-free undamped natural frequency of 1th elastic node 

(rad/sec) 

1th elastic mode normalized mode shape 

a^4(x,y) 

slope of 4>^(x,y) with respect to x (1/ft) 



1th elastic mode motion-dependent generalized force In z 
direction (ft/sec^) 



1th elastic mode gust-induced generalized force In z-directlon 
(ft/sec2) 


INTRODUCTION 

Recent work with control -configured vehicles (CCV) end active 
control technology (ACT) has Improved the performance, stability, and 
handling qualities of large flexible airplanes and has opened up a new 
realm of design frontiers*. 
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With increased sl/e of present day airplanes, and with the increased 
utilization of lighter materials, the elastic behavior of these vehicles 
is becoming an appreciable Influence In their handling and ride qualities. 
Due to the potential adverse effects of elastic mode Interaction with the 
rigid-body dynamics, there Is a need for a simplified method of modeling 
the dynamic aeroelastic equations of motion for use in preliminary 
control system design stages of new airplanes. 

Usually, only calculated values of the rigid-body aerodynamic stabil- 
ity derivatives are available for the preliminary design from sources 
such as DA1C0M2, and little, if any, Infonnation on the stability deriva- 
tives due to elastic modes is available then. However, calculated values 
of the symmetric and antisymmetric orthogonal elastic vibration mode 
shapes and natural frequencies are usually available at the preliminary 
design stage for use In equations of motion formulation. 

We have developed a unique formulation of the equations of motion 
for elastic airplanes that makes use of rigid-body aerodynamic stability 
derivatives and the elastic mode shapes and frequencies to describe the 
aerodynamic forces and moments due to the elastic motion of the aircraft. 
There '*s no need for unsteady aerodynamic theories or experimental data 
on elastic mode aerodynamics as with conventional formulations. 

This paper describes the longitudinal dynamic formulation and veri- 
fication of its accuracy using the B-1 aircraft dynamics at a high sub- 
sonic flight condition. 

EQUATIO NS OF MOTION FO RM! ' LATIO N 

Since the elastic modes do not produce significant drag force 
perturbations compared to the rigid-body motion, the longitudinal 
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equation of motion Is omitted and only the plunge and pitch rigid-body 
equations (short-period approximation) are Included In what follows. The 
formulation of the snail perturbation aerodynamic forces and moments Is 
based on the local effective angle of attack, u(x,y,t), or effective 
plunge velocity, w(x,y,t). where w(x,y,t) • Uon(x,y,t), and the local 
effective pitch rate, e(x,y,t). 

The elastic vibration characteristics are based on the usual approach 


of Idealizing the structure to a flat plate In the xy-plane, and the 
symmetric orthogonal free-free elastic vibration mode shapes <i^(x,y) are 
functions of x and y coordinates In the x,y,z body axes system located at 
the center of gravity^. The sign convention for the mode shapes, mode 
slopes, and generalized displacements Is given In Figures 1 and 2. 

The two time domain short-period and n elastic mode small perturba- 
tion equations of motion about a trim condition are given by (1). 

w(t) - UoO(t) « J w(x,y,t)dxdy + 6^(t) + 

y X ® 


w(x,y,t)dxdy + 


wy w(x,y,t)dxdy 



(1*1 ,2,. . . ,n) 


Zw* , Mq, M^, and are the rigid-body dimensional aero- 
dynamic stability derivatives defined in (2). Z and Z- are assumed to 

q w 

be negligible and are not included In (1). 
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(2) 


The Integral terms and the generalized force terms 1n (1) are 
functions of w(x,y,t), and e(x,y,t), which can be closely approximated by 
(3) and (4). 


n , n , 

w(x,y,t) » w(t) + I ^.{x.v)L(t) - I Uo*.(x,y)5.(l) 

1-1 ^ ^ 1-1 ^ ^ 


e(x,y,t) - e(t) - [ ♦!(x,y)r.(t) 

1-1 ’ ’ 


(3) 

(4) 


The Integral terms can be written as In (5) _nd (6). 


3^Z n 

y X 1-11 


f 

w(x,y,t)dxdy + J - 55 -jy w(x,y,t)dxdy 


f 


y X 


y X 




3 ^ 6(x.y, t)dxdy - M^w(t) + M^w(t) + M^e(t) 


y X 


n p 

I M ^i(t) 
1-1 l-H ^ 


+ M- 5.(t) + M 
^1 ’ ^ 




( 6 ) 


The expressions for Z , Zi , , M; , and MV are tabulated In Appen- 

^1 '1 '1 ‘•1 *1 

dix A. 
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The expression for the mot Ion -dependent generalized force term in 
the n elastic mode equations of motion of (1) 1$ given by (7). 



1 


j Tx l y '^ ♦^(x.y)dxdy 

y X 


where 

Z(t) • M[2„ w(t) ♦ J [Z C.(t) ♦ Z- L(t)] ♦ Z a (t)] 

J-1 J J ® 

and is the ith mode generalized moss. 

Putting (8) into (7), 


(7) 

( 8 ) 



F. , F. , F. , F^ , and F. are tabulated in Appendix A. 

h 'e ^ 

The generalized force term due to C.G. referenced vertical gust 
velocity Wg(t) Is given by (10). 



M 


J 3x3y *i(' 


y X 


Wg(t) 


^ «g<‘> 
“9 


(10) 


Substituting (5). (6). (9), and (10) into (1), Laplace transforming, 
and putting in matrix form yields (11), where four elastic modes (n»4) 
have been explicitly Included. 


3 



9 


STABILITY DERIVATI VEJ 

Since elastic mode shape and slope data Is usually ^Ivun as a func> 
tion of lumped mass stations, which themselves are given in xy-coord1» 
nates, the double Integrals In the tent)s of Appendix A can be conveniently 
represent i as summations over Incremental areas (&x4y) associated with 
each lumped mass point. Thus, It Is necessary to develop methods for 
evaluating the following partial derivative terms at each point: 


!!!!•!!!!'!!!!- !!!!3 

3x3y ’ 3xay ’ 3xay ’ ax3y ’ 3x3y 


Using (2), these become (12) through (16). 

32c. 


32Z, 


8x3y 

3x3y 


-pUqS 


3X3y 


lx3y ' 


?!fw 

-0U0S 



*3x3y 


3x3y 

*3x3y 

S2M 

w 

^2r 

pUoSc \ 


Dx3y 

2t^ 3x3y 



a?r 

pSc 2 m- 


41^ 3x3y 


pUqSc 

TC 


2 3^C, 


m. 


■ix3y 


( 12 ) 

(13) 

(14) 

05) 

(16) 


C, , C„, C , C . . and C, are the total -airplane rigid-body non- 

o a a g 0^ 

dimensional stability derivatives, which are known constants for a trim 
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flight condition. We need to determine the xy area distribution of these; 
I.e., the second partlals In (12) through (16). 

For convent1ona1>ta11*d airplanes, the lift curve slope can be 
reasonably approximated by 


'L " 
a oy 


(17) 


where C. and C, are the wing and horizontal tall contributions. 

‘•ay ‘•Uyj 

Fuselage lift Is neglected as small. Methods for computing these can be 
found In reference 2. The tall contribution Is about ten percent of the 
total. Thus, 


- 0.9 C 


L 


•w 


(18) 


0.1 C. 


‘HT 


A crude approximation, but nne found to be adequate for this forniulatlon, 
Is to assume the derivatives to be uniformly distributed over the xy- 
plane representation of each component (I.e., wing, tall, fuselage). 

More accurate elliptical lift distributions were tried, but resulted In 
very little difference to transfer function dynamics obtained from (11) 
over that for the uniform distributions. Thus, 

r-r* “ 0, for fuselage stations x (y*0) 

(jX oy 




- C| , for wing stations x,y 
w ‘’a 


( 19 ) 


“HT 

^HT 


C, , for tall stations x,y 
^HT '■a 


n 


whtre and art wing and horizontal tall x.y planform areas. 

Cp << and, therefore, Cp will be neglected In evaluation of (13). 

Since C. Is due only to elevator deflection, 
e 


J2C 


L 

Tiiiy 


- • ~ . for tall stations x.y 

^HT 


( 20 ) 


•> 0 , for wing and fuselage stations x,y 


Neglecting small effects due to the fuselage. 


; • C ♦ C 

'm m m 

<1 ay 


a^c 

m 

0 

ax ay” 


m. 


, for wing stations x,y 


m 


• — , for tall stations x.y 


■ 0 , for fuselage stations x (y"C) 


For the pitch damping derivatives, 


C ■ C ♦ C 

m* m* m* 

a Oy 


C - C + C 
^ ^WF '^MT 


(21) 


( 22 ) 


(23) 

(24) 


where WF Indicates wing and fuselage combination. 
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m* 


5xi>y 
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for wing sutlons x.y 


for t«11 stotlons x,y 



3x»y 


m 


‘^WF 


fn. 


M 

‘HT 


for fuselage stations x ,y*0) 


, for wing and fuselage stations x,y 


for tall stations x,y 


(25) 


(26) 


Methods for est1mat1?ig , C„ , C., , and C_ are given In 

’ m* nj* ni„ 


■’HT 


"WF 


'*HT 


references 3 and 4. 

Knowing the airplane elastic mode shapes, slopes, and five rigid- 
body total-airplane stability derivatives, all the terms In Appendix A 
and thus the coefficients In (11) can be computed. 


VERIFICATION WITH B-1 AIRPLANF DYNAMICS 

To verify the accuracy of the unique formulation of elastic airplane 
small perturbation dynamic equations of motion developed above, the tenris 
In equation (11) arc calculated by this method fur the B-1 at a sea level, 
Mach 0.85 flight condition and compared with the corresponding terms In 
equations provided by Rockwell In references 5 and 6, which were generated 
by other methods. 

The xy-pldnform and incremental area divisions for each mass point 
for the B-1 is depicted in Figure 3. Since the ela tic modf*s used In 
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longitudiral dynamic equations are symmetric about the x>ax 1 s, only half 
of the planfo*^ Is shown In the figure. 

As an example of how ^lie double Integrals are evaluated, consider 
the term for the generalised force of the second elastic mode due to the 
third ffode. From Appendix A, It Is 

^ n ♦-(x.y)<;(x.y).ixdy (27) 

y *> 

One half of the value Is obtained by summing ever the 97 stations In 
Figure 3; the other half coming from the symmetric right side planform 
not shown. 

II M 3^2 

The '*erm Is the area associated with each lumped mass point in 

Figure 3. > 2(0 And ^ 3 ( 1 ) are the values of the second mode shape and 

third mode slope In the chordwise x-directlon at the 1 th mass point. 

has three constant values; one for fuselage stations, one for wing 

t»xjy 1 

stations, and one for horizontal tall stations. « 3.94 for this B-1 

flight condition. Wing and tall areas are S ^ * 1946 ft^, • 502 

ft^. From (19), 


0 , fuselage stations 

18.42 x 10“''/ft^, wing stations (29) 

7.14 X 10”‘^/ft^, tail stations 

From (13), with neglected as small. 




ax ay 


a < 
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0, fuselage stations 

-5.708 X 10”'*/ft^-sec, wing stations (30) 

-2.213 X 10'‘*/ft2-sec, tall stations 

The calculation of (28) gives r 2 . ■ 5.8459, which compares with 6.6257 

t.5 

obtained fi*om Rockwell's formulation of the equations. 

All other coefficients In (11) were similarly evaluated and are 
tabulated, along with the values from Rockwell's B-1 equations of motion, 
in Appendix 8. Approximately 80 percent of the terms show good agreement 
with the B-1 data. !n view of the approximate nature of the formulation, 
this Is reasonable confirmation of the validity of the method. 

A further check was made by comparing the roots of the character- 
istic equations for the 13-1 data and this formulation by expanding the 
determinant of the 6x6 matrix of polynomials and coefficients In (11). 

The coupled frequencies In rad/s and damping ratios were calculated for 
each pair of complex roots. The comparisons are shown In Table 1. It Is 


'dxay't 


TABLE 1 FREQUENCIES AND DAMPING RATIOS 
B-1 Data Unique Formulation 


Frequency 

Damping 

Frequency 

Damping 

2.868 

0.489 

3.103 

0.492 

13.298 

0.053 

13.709 

0.034 

21.375 

0.031 

21.221 

0.025 

22.020 

0.020 

22.020 

0.020 

22.480 

0.206 

25.366 

0.233 


evident from the data In this table that the new formulation of the 
equations of nratlon Is surprisingly accurate considering the level of 






IS 


approximations niade. The four symmetric elastic modes of the had 
free-free undamped natural frequencies of 13.591, 14.123, 21.198, and 
22.056 rad/s. All had 0.02 structural damping ratios. The first line of 
numbers In Table 1 corresponds to the short-period frequency and damping 
ratio. 

CONCLU SIONS 

The unique method of formulation of the longitudinal small perturba- 
tion equations of motion for elastic airplanes described herein allows 
the expression of aerodynamic forces and moments, due to elastic vibra- 
tion, In terms of rigid-body aerodynamic stability derivatives. Thus, It 
serves as a useful preliminary design tool for airplane stability and 
control, handling and ride qualities, and control system design studies. 

The good accuracy of the method has been established by comparison 
with more accurate data for the B-1 airplane. The lack of complete 
Information on tne planform geometry of the B-1 and our having to analyt- 
ically calculate mode slopes by curve fits from the mode shape data, 
probably accounts for much of what differences appear In the term by term 
comparisons In Appendix B. Therefore, the new method Is probably even 
more accurate than this one erample comparison Indicates. 

W» have developed a similar formulation for the lateral -directional 
dynamics with the antisymmetric elastic modes Included and are presently 
checking Its accuracy. This formulation will be published In a sequel 
paper to the present one. 
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A PPENDIX A. EQUATIO N COEF FICIENTS 
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/IPPENOIX 3. 

COEFFICIEM VALUES 


Term 

From B-1 Equations 

From Unique Formulation 

''■w 

-0.77480 

-0.69335 


1.3590 

1.4180 


0.80586 

0.81559 


1.7902 X 10’’ 

1.8274 X 10"^ 

ISI 

-0.20177 

-0.01798 


2.4702 

1 . 9202 


0.14486 

0.15373 


-4.7412 X 10"^ 

0.11254 


-8.4911 

-6.7715 


90.322 

103.32 


4.3792 

-1.3083 


-4.1323 

-5.3236 

”c> 

0 

-0.72033 X lO'** 

”C3 

0 

-3.6129 X lO'** 


0 

1.4315 X lO’** 

”c. 

0 

0.06947 X 10“* 

"1. 

-7.5404 X 10“^ 

-10.823 X 10"^ 


50.866 X 10"^ 

2.9502 X 10"^ 

”^3 

7.0361 X 10”’ 

16.427 X 10"^ 

"c. 

-2.0060 X 10"^ 

-3.0665 X 10“^ 


k!1»roducibility of the 
ORIGINAL PAGE IS POOR 


Term 

From B-1 Equations 

Frotti Unique Formulation 


-0.18b49 

-0.04194 


-0.07903 

-0.01592 


0.20945 

0.07459 


-0.09381 

-0.01065 


-0.06630 

-1.1180 

^‘C2 

-0.28286 

17.818 


0.91557 

-0.07921 

f>u 

-0.19759 

-0.69559 


7.2681 

-31.795 

^‘C2 

-15.897 

643.69 


61.886 

9.640 

'■u 

20.894 

-18.803 


0.23360 

0.34121 


-8.2949 

-10.334 


0.11424 

0.00560 


0.11188 

0.11583 


14.002 

18.740 


-306.000 

-415.280 


6.6257 

5.8459 


11.211 

14.696 


-0.12060 

-0.010897 

^’i2 

3.7604 

0.04020 

^n, 

-0.42578 

-0.18203 


-0.25330 

-0.11246 
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Term 

From B-1 Equations 

From Unique Formulation 


7.0455 

2.6015 


33.993 

-15.509 

^.3 

-7.9516 

-1.5492 


3.4837 

2.0596 


-3.2701 X lo"* 

-2.9290 X 10'“ 


24.031 X 10'“ 

25.469 X 10'“ 


-2.7776 X 10'“ 

-3.4422 X 10'“ 

f“U 

-4,0417 X 10'“ 

-4.2249 X 10’“ 


1.6305 X 10'^ 

-0.04974 X 10’^ 


-9.7878 X 10'^ 

3.6896 X 10’^ 


0.70767 X 10'^ 

0.24563 X 10'^ 


1.3340 X 10'^ 

0.18614 X 10’^ 

f>._ 

-22.296 X 102 

-15.978 X 102 

e 


-2.1741 X 102 

-1.5347 X 102 

e 

^33 

6.1537 X 1q2 

4.3685 X 1Q2 

e 

f-s 

0.11048 

0.06489 

e 

F, 

-0.77350 

-0.69335 


P? 

'w 

1.3567 

1.4180 

g 


0.80450 

0.81559 

g 

'•‘“w 

1.7872 X 10'^ 

1.8274 X 10’’ 

g 
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Fl gurt Captions 

Figurt 1. Fuselage Vertical Bending Sign Convention 
Figure 2. Wing Deflection Sign Convention 
Figure 3. B-1 Hass and Area Distribution 









